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INTRODUCTION 

This  paper  deals  with  asymptotic  properties  of  various  models  of 
age-dependent  branching  processes,  relying  heavily  on  Harris  [-5]. 

Part  I  considers  cell  growth  in  which  a  cell  proceeds  in  a  sequen- 
tieO.  manner  through  n  independent  states,  state  R  with  its  life  die 
tribution  R=l,2,  ...,n.  At  the  end  of  mitosis,  the  state,  the 

cell  divides  into  similar  cells,  the  number  of  -vdiich  is  governed  by  a 
generating  function  h,  independent  of  the  time  and  other  cells  of  the 
process . 

Let  Z(t)  be  the  number  of  cells  at  t,  and  Z^(t),  R=l,2, ...,n, 
be  the  nvimber  of  cells  in  state  R  at  t.  For  increasing  cell  popula¬ 
tions,  that  is,  h*  (l)  s  m  >  1,  we  treat  the  fo3J.owing  topics,  (t^w 
throu^out) . 

(1)  E[Z(t)]  and  Var[Z(t)] 

(2)  forward,  backward,  and  total  times;  Given  that  a  cell  is  in 
the  R^^  state  at  t,  to  find 

(a)  the  distribution  of  its  time  to  reach  the  end  of  that  state, 

(b)  the  distribution  of  the  time  since  it  entered  that  state,  and 

(c)  the  total  time  it  will  have  spent  in  that  state,  respectively. 

(3)  asymptotic  fraction  of  cells  in  state  R[s  P^,]  and  the  relation 
between  1^  and  E[Zj^(t)], 


A  variation  of  the  above  model  treats  cell  growth  in  ’which  the  cell 
proceeds  from  state  to  state  according  to  a  general  semi-Markov  process 
until  the  mitotic  state  is  completed,  when  division  into  similar  cells 
in  accord  with  h  occurs. 

We  treat  the  5  state  semi -Markov  case  in  detail  to  obtain 

(1)  the  equivalent  life  distribution  G  of  a  cell, 

(2)  asymptotic  fraction  of  cells  in  a  state, 

(5)  forward,  backward,  and  total  times  in  a  state. 

For  the  general  n-state  semi -Markov  case,  espressions  for  (2)  and 
(5)  above  are  fovmd,  and  it  is  shown  that  as  h*  (l)  =  m>lfl,  that  the 
asymptotic  fraction  of  cells  in  a  state  approaches  a  simple  limiting 
form  which  agrees  with  a  result  of  Saith  [?]• 

For  both  the  sequence  of  states  and  general  semi -Markov  model,  the 
convergence  of 

ZpCt) 

“  ELzJ't)]  *  R=l,2,  ...,n, 

a  non-degenerate  random  variable,  in  quadratic  mean  is  discussed. 

The  correlations 

lim  p(W^(t),  Wg(t))  =  =  1  , 

t-><» 

R,S  =  1, 2,  ...,n,  and  W  =  W  a.e.  Moments  of  Ws 

K  b 

obtained . 

Bart  II  considers  the  random  variable  R(t)  =  the  total  ntimber  of 
births  by  t  in  a  simple  age  dependent  branching  process  for  m  >  1. 

The  following  topics  are  treated  for  t  large: 


W_  may  be  easily 


(l)  E[N(t)]  and  Var[N(t)3 

a  non-degenerate  random  variable,  in 

quadratic  mean. 

(5)  Let  Z(t)  =  number  of  cells  at  t. 

From  t  5]^ 


Wj^(t) 


in  quadratic  mean  . 


It  is  shown  that  the  correlations  satisfy 


lim  :p(W  (t),  U(t))  =  ;p(w  ,  wj  =  1 

■t-»  00  o  X  ox 

and  that  ¥  =  W,  a.e. 
o  X 

(4)  The  results  for  E[N(t)]  are  checked  with  the  corresponding 
discrete  case  by  matrix  methods. 

(5)  For  the  case  m  <  1,  W(t)  a.e.  Moments  of  may  be 
calculated. 

(6)  For  m=l  and  m  <  1,  moments  of  E(t)  are  discussed  for  both 
discrete  and  continuous  time. 

Bart  III  treats  the  case  of  two  types  of  cells,  of  ■vdiich  only  one 
type  may  divide  while  the  other  either  accumulates  or  is  eventually 
absorbed  in  the  medium.  An  example  of  this  is  the  production  of  stem 
cells  and  red  blood  cells  from  parent  stem  cells.  Two  related  models 
are  considered,  both  examples  of  the  “reducible  case".  For  an  outline 
of  the  irreducible  case,  see  Snow  [8].  Increasing  cell  populations 
are  considered. 
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A  binary  fission  case  of  each  model  may  be  represented  schematically 
to  Indicate  the  various  combinations  of  cell  births  with  corresponding 
probabilities,  where  type  2  represents  the  proliferating,  and  type  1 
the  non-proliferating  cell. 

Model  I 

A 

2^2 
P  2pq  d 

where  p  +  q  =  1, 

Model  II 

P  q 

where  p  +  q  =  1. 

Let  =  number  of  cells  of  type  i  at  t,  i=l,2. 

N^(t)  =  total  number  of  bliijhs  of  type  i  by  t. 

The  following  topics  are  treated  for  t  large: 

(1)  means  and  second  moments  of  Z^(t),  N^(t),  1=1,2. 

(2)  Results  for  the  means  in  (l)  are  compared  with  the  discrete  case 
by  matrix  methods,  and  a  comparison  is  made  with  results  of  Snow  [8  ] . 


k 


(5) 


z.Ct) 

W.  (t)  =  '  -»W.  In  quadratic  mean,  i=l,2 

^  ELZ^(t) ]  ^ 

,  ,  N  (t) 

W  .  (t  J  s  — ' "  -»  W  .  In  quadratic  mean,  i=l,  2, 


(4)  The  limits  of  the  pairwise  correlations,  as  t  -»«,  of  W^(t), 
Wg(t),  W^3^(t),  respectively  equal  to  the  corresponding 

pairwise  correlations  among  W^,  Wg,  W^g,  and  are  all  equal  to 

1,  and  =  Wg  =  =  W^g  a.e. 


5 


PART  I  -  MULTI-STATE  CEIiS 


1.1  Sequence  of  States. 

A  cell  exhibits  growth  and  division  into  other  cells  according  to 
the  scheme  to  be  described. 

The  cell  enters  state  1,  for  which  it  remains  a  time  which  is  a 
random  variable  with  distribution  function  then  proceeds  to 

state  2,  with  distribution  function  F^,  and  so  on  until  it  completes 
the  n^^  state  of  growth,  the  mitotic  state,  after  which  it  divides. 

The  overall  life  distribution  of  the  cell  is 


G  =  F^  *  Fg  . 

The  cell  divides  into  k  other  cells  with  probability  Pj^, 
k=0,l,2, ...  each  with  independent,  identically  distributed  growth 
pattern  as  the  parent  cell,  and  independent  of  each  other. 

Define  h(s)  =  T  p,  s'*^  . 

k=0  ^ 


Schematically: 

H — I - f— 

X, 

A  brief  treatment  of  the  two- state  case  is  given  in  Ch.  6, 

sect.  26,  where  a  formula  for  the  limiting  mitotic  index,  which  is  the 
limiting  fraction  of  cells  in  mitosis,  is  given. 

As  an  example,  certain  normally  proliferating  cells  pass  throu^ 
four  states  to  division,  called  the  G^,  S,  Gg,  and  mitotic  state  in 
that  order.  Each  corresponds  to  a  certain  condition  of  the  cell  and  has 
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a  certain  sensitivity  to  radiation.  For  example,  the  S  state  is  the 
DNA  synthesizing  state,  and  it  along  with  the  mitotic  state  is  more  sen¬ 
sitive  to  radiation  than  either  the  or  Gg  which  are  called  rest¬ 

ing  states.  The  results  of  this  paper,  especially  part  I  may  aid  In 
the  computation  of  the  fraction  of  cells  heavily  damaged  hy  radiation. 
See  Stohlman  [9]. 

1.2  E[Z(t)L  Var[Z(t)]. 

To  discuss  cell  behavior,  it  is  convenient  to  Introduce  the  gen¬ 
erating  function  of  cell  life. 

Definition.  Let  Z(t)  denote  the  number  of  cells  at  time  t, 

starting  with  one  cell  at  t  =  0. 

Definition.  F(s,t)  =  |  P[Z(t)  =  k]s^  . 

k=0 

We  have  [5]  (Ch.  6,  sect.  7)  the  relation 

F(B,t)  =  s(l-G(t))  +  f  h(F(B,t-u))  dG(u) 

-'o 

m(t)  =  E[Z(t)]  = 

m(t)  =  1-G{t)  +  m  f  m(t-u)  dG(u)  . 

-^0 

We  assume  m  s  h’ (l)  >  1,  h"(l)  <  w.  Define  a  >  0  by  the  equation 
e”^^dG(u)  =  By  lenmaa  1  of  the  appendix,  as  t  -^oo. 


m(t) 


[1-G(u)  ]  du 


ue"^^  dG(u) 
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or 


m(t)  - 


m-1 


am 


C  eo 

-0£U 

ue 

Jo 


dG(u) 


e  =  n^e 


A  further  resvilt  has  ( [j]  Ch.  6,  sect.  19)  for  t 


Z(t)  ~  W  n^e 


ext 


where  W  Is  a  random  variable  such  that 


a^(w)  = 


(m+h"(l))  /  e“^‘^^dG(u)-l 

Jo 


r 

Jo 


1-m  /  dG(u) 


To  find  Var[Z(t)]^  note 


Var[z(t)] 


Ss' 


S=1 


and  using  theorem  18.1  of  [3]j  Ch.  6,  we  obtain,  as  mi(^t  be  expected, 


Var  Z(t)  «  a^(V)  n^e^^* 


1.5  Forward,  backward  and  total  times. 

1.5a  Total  cell. 

The  total  backward  time  is  the  distribution  of  time  the  cell  has 
been  living  until  t,  given  that  it  is  alive  at  the  present  time  t. 
Let  A(x)  =  distribution  function  of  the  backward  time.  Then  (tj], 
Cli.  6,  sect.  24) 
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[l-G(t)]  dt 


A(x) 


r  00 

I  [l-G(t)]  dt 

Jo 


We  may  find  the  forward  time  distribution,  denoted  by  B(x),  for 
the  total  cell  by  using  A(x). 


P„(t,x)  =  P  [cell  will  die  by  age  t+x|age  is  now  t] 


l-G(t) 


We  then  have 


B(x) 


Jo 


e"“^  [G(t+x)  -  G(t)]dt 


Pg(t,x)  dA(t)  =  — 


f  e"^  [l-G(t)  ]dt 

Jo 


Denote  the  total  time  distribution  for  the  entire  cell  life  by  C (x) 


P_(t,x)  =  P[cell  will  die  by  age  xlage  is  now  t]  . 


P^(t,x)  =  0  ,  t  >  X  . 


C(x) 


x)  dA(t)  = 


e"^  [G(x)  -  G(t)]dt 


f  e'^*  [l-G(t)  ]dt 

Jo 


Note  A(x)  *  B(x)  ^  C(x). 
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1.5b  Each  state. 

Given  that  the  cell  proceeds  sequentially  through  n  states  before 
mitosis,  with  corresponding  distribution  functions  Fg, ...,P^. 

As  before,  G  =  *  Fg  *...*  F^  and  a  is  such  that 


e-«^  dG(u)  =  i  . 


To  find  A^(x),  the  backward  time  in  state  1,  we  proceed  as  in 
Theorem  2l^.l  of  [5]^  Ch.  6. 

Denote  Z^(t)  =  the  number  of  cells  at  t  in  state  1. 

Z^(t,x)  =  the  number  of  cells  at  t  in  state  1  whose  age 

is  <  X.  Then  if  F.  (s,t)  =  ^  P[Z  (t)  =  jjs*^ 

J=0 


F^  (s,x,t)  =  J  P[Z  (t,x)  =  j]  s‘ 
J=0 


we  obtain 


F^(s,t)  =  [1-F^(t)]s  +  F^^Ct)  -  GCt)  +J^  h(F^(s,t-u))  dG(u) 

E[Z^(t)  ]  s  m^(t)  =  l-F^(t)  +  m  J"  m^Ct-u)  dG(u) 
F^(s,x,t)  =  [l-F^(t)](sJ(x-t)  +  l-j(x-t)]  +  -  G(t) 

+  J  h(F^(s,x,t-u))  dG(u) 


where  J(r)  = 


1,  r  >  0 
0,  r  <  0 


E[Zj^(t,x)]  =  ra^(t,x)  =  [l-F^(t)]  J(x-t)  +  m  J  m^(x,t-u)dG{u) . 
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Solving  for  m^{t,x)  and  in^(t)  for  t  large  by  lemma  1  of 
the  appendix^ 


- 


“  — -  ♦ 

e“^  [1-F^(t)]dt 


To  find  Bjl(x),  the  forward  time  in  state  1^  ve  use  A^(x)  as 
done  before. 

Let  P^(t,x)  =  P[cell  will  leave  state  1  by  age  t+x[age  in  state 
1  is  t  ] 

P  (t+x)  -  (t) 

P^(t,x)  =  * 


\(x)  =  P^(t,x)  dA^(t)  = 


JJ  [F^(t+x)  -  F^(t)]dt 
J  e’°*  [1-F^(t)]dt 


To 


obtain  C^(x),  the  "total  time"  in  state  1,  define 
I-IC  (t,x)  =  P[cell  will  leave  state  1  by  age  x|age  in  state  1  is  t]. 


P^c('tjx)  =  - r-Fj^'(t') - -  '  ^  ^  ^  '  X  <  t 


.00  f  [F^(x)  -  F^(t)]  dt 

=  L  ^1^^^  -vv  ■  - 

®  j  e"^  [1-P^(t)]dt 


Let  A2(x)  be  the  backward  time  in  state  2. 
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Denote  lay  Z2(t,x)the  number  of  cells  in  state  2  at  time  t  -whose 


age  is  <  X. 
Since 


and 


defining 


P[X^  +  X2>t}  t>X^>  t-x]  =  J  /■  dF^Cxg) 


+  Xg  >  t;  t  >  X^]  =  F^(t)  -  *  FgCt) 


Fp(s,t)  =  2  P[Z  (t)  =  Jls*^ 

F  (s,x,t)  =  I  P[Z  (t,x)  =  Jjs*^ 


J=0 


we  obtain 


FgCs^t)  =  s[P^Ct)-F^  *  FgCt)]  +  l-G(t)-[F^(t)-F^  »  F^Ct)] 


/. 


+  /  h(Fo(s,t-u))  dG(u) 


and 


F2(s,x,t)  =  s[[^  dF(x^)[  dF2(x2)]  +  l-G(t)-[f  ^  dF  (x^)f  dFg(xg)] 

-'t-x  t-x^  -^t-x  "^Jt-X^ 

L 


+  1  h(F„(6,x,t-u))  dG(u) 


dF^Cs^t) 


nigCt)  s  EEZgCt)  ]  - - ^ 

,  %  - 
TUgCt^x)  =  EEZgU^x)]  = 


S=1 


3Fg(s,x,t) 


s=l 
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Using  lemma  1  of  the  appendix,  for  t  -»oo,  we  find  that 


m2(t,x) 

in2(t) 


~  AgCx)  = 


*  Fg(t)]at 


To  obtain  BgCx),  the  forward  time  in  state  2,  define 


^  ®2U2)lat 


P2^(t,x)  =  P[t+x  >X^  +  X2>t  i  X^<t] 
rt  rt+x-X 

PgjCt,*)  =  ’■dFgCx^). 


Then 


=  /“ 
-*0 


BgU)  =  P^^(t,x)6A^(t)  = 


r'^  f  t+x-x. 


e"°*[F^(t)-F^  *  F2(t)]dt 
Similarly,  for  CgCx),  the  total  time  in  state  2,  define 

P2^,(t,x)  =  P[x  >  X^  +  ^2  >  t  j  X^  <  t],  X  >  t  i  P2f,(t,x)  =0,  t  >  X. 


r  t 


i"  x-x. 


P2g(t,x)  =  dF2(x2)  ,  X  >  t  ;  p2j^(t,x)  =  0  ,  t  >  x 


Jt- 


^gCx)  —  J*  P2Q(t,x)dA2(t)  = 


e"°^tF^(t)-Fi  *  F2(t)3dt 


15 


In  general,  to  find  Aj^(x),  Bj^{x),  and  Cj^(x),  k  >  2,  define 
C^(x)  =  *  Fg  *•••*  Fj^(x),  then  we  use  the  same  arguments  as  for 

A2(x),  B2(x),  C2(x),  where  corresponds  to  F^^  and  F^^  to  Fg 

Thus  we  obtain,  for  k  >  2 


-  ± 


\M  = 


dFj^(xg)jdt 


Jq  ®  *^[^_lC't)-\(t)]dt 


C.  (x)  = 
k'  '  roo 


/• 

•^0 


It  is  easily  shown  that  all  of  the  above  distributions  may  be 
derived  solely  by  means  of  generating  functions  following  theorem  214-.  1 
of  (t3],  Ch.  6),  as  was  done  above  in  obtaining  A^(x)  and  AgCx). 

Further,  it  is  again  possible  to  derive  all  of  the  above  distri¬ 
butions  very  guickly  by  heuristic  arguments  similar  to  one  due  to 
R.A.  Fisher.  See  ([5],  Ch.  6,  sect.  2l|-). 

1-^  AsyTtg)totic  fraction  of  cells  in  each  state. 

As  before,  we  consider  cell  growth  through  a  sequence  of  states  as 
defined  above,  making  the  same  assumptions  and  using  the  same  notation. 


)]dt 
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Definition;  Pp(y)  =  P[cell  of  age  y  Is  in  state  li],  R=l,2, 
Definition;  P_  =  asymptotic  fraction  of  cells  In  state  R. 

A 

For  the  case  n  =  2,  P^  is  obtained  in  [3],  Ch.  6,  sect.  26.  In 
general,  P^  is  called  the  mitotic  index.  ¥e  have 


PrCf)  =  P[  \  <  y; 


1=1 


n 


I  X  >  y] 


1=1 


As  before,  denoting  =  F^  *  Fg  *•••*  P^(x)  and 

G(x)  =  F^  *  Fg  F^  we  may  write 


Pj,(y) 


^g0?1 


*■5  =  _(  \(yWy)  -  If 

1.5  E[Zp(t)]  ,  Var[Zj^(t)]. 

We  now  obtain  the  asymptotic  mean  and  variance  of  the  number  of 
cells  in  each  state. 

Let  =  number  of  cells  in  state  R  at  t . 

Define  Pj^(B,t)  =  2^  P[Zj^(t)  =  Jjs'^. 

Dy  the  arguments  in  the  previous  section,  we  obtain 

Fj^(s,t)  =  s[Gjj_^(t)-Gj^(t)]  +  l-G(t)-[Gg_^(t)-Gg(t)] 

+  J  h(Fj^(s,t-u))dG(u)  . 
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mj^(t)  s  E[Zg(t)]  =  — ^ - 


s  B=l'  using  lemma  1  of  appendix 


2  for  t  large. 


m^(t)  ~  — - - 

m/  te'^dOCt) 


e®  =  nj,e°* 


3  F  (s,t)| 

Also,  since  Ver[Zj,(t)  1  =  — S-g—  .[yt)]^for  t  lerge, 

OS  8=1 

using  the  method  of  lemma  18.1  of  [3],  Ch.  6,  ve  obtain 


Var[Z^(t)]  ~  r? 


(h»(l)+m)  fV^^GCu)-! 

_ £o _ _ 

1.^  r  dG(u) 

JO 


As  noted  in  [3],  Ch.  6,  sect.  l8,  the  denominator  is  positive  since 


[*6**^“^  dGCt)  <  r, 

-•'o  Jo 

Prom  the  first  section  of  this  paper, 


e""*"  dG{t)  =  -i  . 

m 


E[Z(t)]  s  mCt)  = 


o  r“ 

am  /  ue‘^dG(u) 

Jo 


Hence,  as  we  may  have  expected,  for  t  large, 

mjj(t) 
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Hie  aTaove  quantities  are  of  little  practical  interest  in  the  cases 
m=l  and  m  <  1,  since  it  is  known  that  the  populations  die  out  with 
probability  one  in  those  cases.  See  [5],  Ch.  6,  sects.  1-lk. 

1.6  Cell  growth  by  states  according  to  a  semi -Markov  process. 

We  consider  the  following  model.  A  cell  evolves  via  states  to 
mitosis,  but  instead  of  proceeding  sequentially  from  one  Intermediate 
state  to  the  next  in  a  deterministic  way,  instead  the  state  selection 
proceeds  according  to  a  Markov  chain.  Specifically,  given  n  states. 
Including  the  mitotic  state  (denoted  by  the  n^^  state),  assvime  there 
exists  an  n  x  n  irreducible  positive  recurrent  transition  probability 
matrix  P,  with  zero  trace,  giving  the  probabilities  of  transition  from 
one  state  to  the  next  during  the  growth  of  a  cell.  As  soon  as  the  n"^^, 
or  mitotic  state,  is  completed,  the  cell  divides  into  r  cells  with 
probability  Pj.,  r=0,l,2, ...  each  with  a  growth  pattern  independent 
and  identically  distributed  as  the  parent  cell,  and  Independent  of 
each  other .  Define 

00 

h(s)  =  3  P  s^  ,  h*  (1)  s  a  >  1  . 
r=0  ^ 

The  time  spent  in  the  k^^  state,  k=l,2,  ...,n  given  that  the 
next  transition  is  into  state  J,  for  ,J=1,2, ,..,n  and  k  /  j,  is 
a  random  variable  with  non- lattice  distribution  function 

Fj^(O)  =  0,  dependent  on  the  _k^^  and  states,  but  otherwise  inde¬ 
pendent  of  the  state  of  the  system. 

At  present,  thereappears  to  be  no  physical  example  of  this  process 
in  cell  growth,  but  a  possible  Interpretation  is  that  a  cell  which 
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evolves  from  state  to  state  selects  the  most  "accesslhle"  state  at  each 
transition,  vhere  "accessihility"  may  be  considered  to  be  determined  in 
accord  with  a  Markov  process. 

We  begin  with  a  study  of  the  3-Btate  semi-Markov  model.  For  this 
case  we  compute  the  equivalent  life  distribution  G(t)  of  the  cell, 
the  forward,  backward,  and  total  time  distributions  for  a  particular 
state,  and  the  asymptotic  probability  of  a  cell  being  in  a  given  state 
[asymptotic  fraction  of  the  culture  in  a  given  state].  Further,  we  let 
mil,  or  equivalently,  a 4,0,  for  a  defined  by 

r  e"‘^%G(u)  =  ^  , 

Jo  "" 

and  compare  our  results  with  those  of  amlth  [7]- 

To  generalize,  for  the  case  of  n  states,  we  compute  the  asympto¬ 
tic  fraction  of  cells  in  a  given  state  as  a4'0,  and  again  compare  with 
Smith  [7 ] . 


1.7  Three  state  semi -Markov  model. 


Let 


Schematically, 


a  b^ 
o  d  , 
q  oj 


a+b  =  c+d  =  p+q  =  1 
a,b,c,d,p,q  ^  0 
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state  "5"  is  the  mitotic  state  in  the  figure.  Define 


=/g  e-‘yaFy(y) 


k=l,2,3 

^=1,2,3  and  J  k. 


l.Ta  Equivalent  life. 

A  consideration  of  the  various  ways  in  which  a  cell  may  progress 
through  the  mitotic  state  yields  the  following  Laplace-Stleltjes  trans¬ 
form  for  the  equivalent  life  dlstrlhution  G. 

OO 

4q  =  (P'b  pda  4y^2.-^  qcb  <ld  q  4y^ 


j  (plJ  pda  ^21^13^  ^^23^  ^5l'^  ^  ^32^ 

1-ac  021^12 

1.7b  Fraction  of  cells  per  statej  Pjj. 

Let 


=  asymptotic  fraction  of  cells  in  state  1 
q^fy)  ~  P[cell  in  state  1  at  age  y]. 

Since  the  equivalent  age  distribution  of  the  cell,  A(x),  is 


A(x)  = 


Pe"'^[l-G(t)]dt 

Jo 


r 

^0 


e"°*[l-G(t)]dt 


we  see  that  ([3]>  Ch.  6,  sect.  2€) 
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where  a  >  0  is  defined  by 


To  compute  q^(y)^  write 

Ql]^(y)  =pPEcell  in  state  1  at  age  yjcell  born  into  state  1] 

+  qP[eell  in  state  1  at  age  y|cell  born  into  state  2], 

so  that 

00 

+  ac  '“^»P2i(y)-(F;^g*P2l)  '“>»P2j^»(aFj^2+bFj^j)  (y)  ] 

where  Is  the  convolution  of  F^^  ,  i=l,2,  and  J=l,2,  3^1, 

and 

=  uU)  = 

We  have 

(ot) 

r«*  ^ 

i  =  -V 

^00 

where  e"^^dq^(y),  the  Laplace-Stieltjes  transform  of 

^l^ty)  evaluated  at  a. 

Ihking  Laplace-Stielt Jes  transforms  evaluated  at  using  the 
previous  expression  for  q^(y), 


1  I  >  0 

« 

0  e  <  0 
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n=0 


12  “=^15' 


+  qo  X  (ac)“4^Wj^»i2l)“tl-(as;;^3+ W„)J  (a) 


n=0 


12  “^^13' 


Since 


L  («) 


a 


l“(asifi2+  (a) 

P  +  q,c  9*'2i(Q!)  ] 

d  J 

1-ac  J 

l-^gCa)  f«> 


a 


=  J”  [l-G(y)]e-°^ay  =  i[l-i] 


p  a 

■  ~Wg(a7 

a 


’l-Cajif^g+b^^^)  (a)' 

p  +  (ic/^2j^(a) 

a 

Ll-ac  ^3^2 

-  [1-  -] 

a 


For  comparison  with  to  be  defined,  we  express  in  terms 


of  1=1, 2, 5,  j=l,2,3,  ±f:  j.  Then 


P,  = 


I"  (a^i2'''^^15^ 


a 


[p  + 


L-ac^^2  ^21  +pdaj^^2(o=)  5^23  ^“^'*‘*^‘^^^21  ^13  ) 


(p^^^(a)+qj^^2(‘^^) 


l-7c  Comparison  with  a  result  of  Smith. 

To  compare  this  with  a  result  by  Smith  [7]y  we  compute 


/: 


[1-  (aP^2  (y)+bF^^  (y)  )  ]e~^^  dy 


^1  fO° 

/  [1-Ki(y)]e"°‘ydy 

Jo 
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where  R^Cy)  is  the  recurrence  time  distribution  of  the  event; 

[cell  branch  enters  state  1].  That  is,  we  follow  any  cell  branch,  or 
strand,  until  state  1  is  re-entered  for  the  first  time  in  order  to 
compute  Rj^.  We  mention  that  an  unsolved  problem  is  to  follow  all  cell 
branches  until  state  1  is  re-entered  for  the  first  time  among  all 
branches,  and  to  compute  a  recurrence  time  distribution  R(y)  for  this 
case. 


where 


=^1=  - - 

a 


r  _ _ _ _  Q=.  _ 

l-(jd9('2^^^2'' (®®i^2_2^2l'*'^^^^12^25^3l"*'  ^13^52^21"**  "^^^13^31^ 

a(l-qd^2^9[^2) 

In  general,  for  a  >  0,  ^-^4  definition  of  recur¬ 

rence  time  distribution  R,  Smith's  theorem  5  [?]  can  not  be  extended 
to  this  semi-Markov  process  in  this  way. 

We  may  let  a  4^0  in  the  expressions  for  and  r^. 

Setting  a  =  0,  we  see  that 
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yielding  Smith's  result  in  the  limit  as  mil.  Puorther,  for  the  case 


a  =  0, 


f 


udG(u) 


vhere 


5  r“ 

kij  =  X  J  and 

ih 

3 

£  =  (rt  ,  jT,)  ,  ^  n  =  1,  rt  >  0  ,  i=l,2f'^f  with  the  relation 

±  C.  0  1  1 


JT  =  It  p,  in  this  case 


(it  ,itp,ir  )  =  (it  ,it  ^it  )  /  O  a  b\  . 

"■'"^IcOd 

\P  qo/ 

It  may  be  shown  by  a  general  theorem  to  follow  In  a  later  sectlonj  that 
for  a  =  0^ 

■P  ^l^^l 

To  complete  the  discussion  of  the  5-state  semi-Markov  model,  we 
give  the  backward,  forward,  and  total  time  distributions  of  state  2. 
There  are  computed  using  the  previously  derived  formulas  for  the  dis¬ 
tributions  in  the  sequence  of  states  model  and  the  law  of  conditional 
probability.  We  sum  up  the  backward,  forward,  and  total  time  distri¬ 
butions  respectively,  each  corresponding  to  a  possible  path  ending  in 
state  2,  with  each  weighed  by  the  asymptotic  probability  or  frequency 
of  that  path,  given  that  the  cell  eMs  up  in  state  2.  Analogous  expres¬ 
sions  hold  for  state  1. 
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1.7d  Forward^  ’backward,  and  total  times. 

Let  P'2  -  asymptotic  fraction  of  cells  In  state  2. 

Define,  for  n=0,l,2, ... 

a^(y)  =  pa(ca)”[F^2*(^12*^21^  (y)-F^2* ^^12*^21^  ^ 

so  that 

a^Cy)  =  P 


“cell  of  age  y  is  in  cell  born  in 
state  2  for  n^th  time  state  1 


and 


f  ri-G(y)  ]e"^^dy 

Jq 


cell  is  in 
state  2  for  n^^  time 


cell  born 
in  state  1 


Define,  for  n=0,l,2, ... 

■fci^(y)  =  aCac)*^!  (F^2*^21^  “  ^^12*^21^  ^ 


so  that 


b 

n 


(y) 


cell  of  age  y  is  in 

cell  born  in' 

^  state  2  for  n^^  time 

state  2 

and 


h  = 
n 


dy 


2k 


=  =4 

n  k 


cell  is  in 


cell  "born 


Then 


^n  ”  ^  state  2  for  n^^  time  in  state  2  * 


^2^1  ^ V^n^  * 

n=0  ^  ^ 


We  then  have 


/  6-“*  [Pj^2*(F3^2*F23^)  (t)-Fj^2*(Fj^2*^Sl)  '“'*(oF23^+aF25)  (t)  Jdt 


Ag(x)=  ^  2  a 
2  n=0. 


r  d[cF2^+dF23](x2)]dt] 


/o  ®''^f^l2*f^21*^12)  (^5-^12*f^l2*^21^  ^^^*(cF2^+dF32)  (t)  Jdt 


r^2*(Fi2*F2l) 


^\"^[(F^2*^2i5  (t)-(F^2*^2^)  ^'^^^(cFg^+dFg^)  (t)  Jdt 


k  nio""  ^ /tl^^  ^  (-l)  f  ^t'=^2l"^255  ^  ("2) 


e"°*  t  (^12*^21^  ■  ^^12*^21^  ^''^^(cFg^+dFg^)  (t )  Jdt 


p  J^^^*fcF  +dF  1 

12  ^21''  ^*^^21  *^25^ 


„  £^"'^it=J'2i-*dF25](x2)]dt 

B2(x)=  ^  2  ®n  "Tcq  - - 
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J  [  (cPg^+dFg^)  (t+x)-  (cPg^+dFg^)  (t)  ]dt 

r  00 

/  e"^[U(t)-(cF2^+dF2^)(t)]dt 


T 

f: 


■2  n=l 


J  d[(F^2*F2^)^^^](x^)J'  ^d[cF2^+dF2^]  (x^)  ]dt 


where 


U(t)  = 


Ij  t  >  0 

0,  t  <  0 


0£(k)=  a  I  a  , 


2  n=0 


Jo^'^ (Xj)  ]« 


Jo 


e"“^[Fj^2*(Fi2*F2i)  (t)-Pi2*fFl2*^21^  ^^^*(cF2i+dF23)  (t)  ]dt 


Jo 


e"*^  [  (cFg^+dFg^)  (x)-  (cFg^+dFg^)  (t)  ]d-t 


[U(t)- (cFg^+dFg^) (t) ]dt 


-at 

e 


+  ^  S\ 

2  n=l  “ 


Jo 


[J^  d[ (F^g^Fg^)  ]  (x^)  j^  %[cFg^+dF2^]  (Xg)  jdt 


r  e"^  [  (Fi2*^21^  ^^^'^^21*^12^  ^^^^(cFg^+dFg^)  (t)  Jdt 

0 


We  return  to  the  discussion  of  the  asymptotic  fraction  of  cells  In 
a  state,  hut  this  time  consider  the  semi-Markov  case  for  an  arbitrary 
number  n  of  states  which  are  available  to  the  cell  in  its  evolution. 


26 


1.8  n-state  seml-Marlov  process;  Pj^Co:)  as  ocJ,l. 

Theorem.  Let  P  "be  an  n  x  n  irreducible  transition  matrix  with 
zero  trace .  Let 


P,  (a)  =  asymptotic  fraction  of  cells  in  state  k,  k=l,2, ...,n, 

A 

where 


G  the  equivalent  cell  life  distribution. 


P- 


k 


n  roo 

J  Pfej  j  k=l,2,  ...,n, 


where 


<  jr,  satisfy  ic  >  0,  k=l,  2,  ...^n,  T 

A  A  A 

j  5  JC  —  itP* 

Then,  as  a-Vo  , 


n  ^  k=l,2,  ...,n. 

1 


s=l 


s  s 


Proof. 


P,  = 

k  z'  w 


qjj(y)e"^dy 


where 


(  [l-G(y)] 

Jo 


e”^dy 


q^(y)  =  P[cell  at  age  y  is  in  state  k], 
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Letting  oc  ^  0 


"K 


f  [l-G(y)  ]dy 
•Jc 


For  Gg  continuous  distribution  functions  on 

G^(0)  =  GgCo)  =  0,  it  is  itnmedlate  that 


^00  ^eo 

*  G2(y)]dy  =J^  udG2(u) 


Hence 


r 

Jo 


=  I  ihr  =  — 

y.jj.v-jf-'uy  ^  n  me  ^  "Jt 
u  r=-L 


n 


vhere 


p(r) 
n  nk 


ra  cell  considered  as  a  Markov  chain  hits 
ij  state  k  at  r"^^  step  without  hitting 
Lmltotic  (ntb)  state 


See  Chung  [2],  pp.  h'^-kh,  h9. 
Hence 


I  [l-G(y)  ]dy  =  2  W 

Jo  n  k=l 


and 


P.  = 


Vk 


k  n 

£=1  ^  ^ 


Tbe  theorem  is  proved. 


[O,  w)  with 


cell  born' 
at 

0'*'^  step  . 
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Two  observations  are  in  order.  The  quantity  is  Interpreted 

n 

as  the  expected  visiting  time  to  state  k  before  mitosis,  the  n't!!  state 


[2]. 


We  write  P,  ,  for  a  =  0,  as 
k 


ii  +  \  ■'  *11  "ii 

^  jTi  ^ 


The  quantity 


n  It  |i 

T  -iJ. 


is  then  the  expected  time  spent  in  states 


J7k  \ 

other  than  k  before  returning  to  k. 

00 

The  theorem  holds  for  a  countable  infinity  of  states,  if  T  n.ii  <<», 

k=0  ^  ^ 

•v^ere  state  0  is  the  mitotic  state.  This  form  for  Pj^  yields 
Smith’s  theorem  5  [7]  Toi*  our  model. 


1.9  Convergence  in  quadratic  mean  and  pairwise  correlation. 

To  conclude  Bart  I,  we  investigate  the  asymptotic  behavior  of  the 
random  variables 

Zj,(t) 

=  E[Zjj(t)  ]  R=1.2,...,n 

for  the  general  seml-Jferkov  model,  which  includes  the  sequence  of 
states  model  as  a  special  case. 

Theorem.  Let  m  >  1,  and  h"(l)  <  m.  In  the  semi -Markov  model, 

Wj^(t)  f  R=l, 2,  ...,n, 

a  non-degenerate  random  variable,  in  quadratic  mean. 

The  pairwise  correlations  satisfy 

Ito^  p(Wp(t),  Wg(t))  =  p(Wj^,Wg)  =  1,  R,S=l,2,...n  . 
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Furthermore, 

Sketch  of  I^oof.  See  Part  II  for  a  aimilar  proof  carried  out  in 
complete  detail.  Major  steps  are  given  here. 


Let 


Gp(s,t)  =  2  P[Zp(t)=j]s'^ 


j=0 


Let 


Po(t)  =  P[cell  in  state  R  at  t] 


Let  G(t)  =  equivalent  life  distribution  of  the  cell. 
Then 

Gj^(s,t)  =  s  Pj^(t)  +  1-G(t)-Pj^(t) 

h(Gj^(s,t-u))dG(u) 

By  previous  methods,  as  t  ->#9,  since 


E[yt)] 


5Gj^(s,t) 

5s 


s=l  ' 


•where  a  satisfies 


and 


50 


Further 


2  S  GpCsj't) 

ElzJCt))  = — s_^ 


as‘ 


S=1 


and  we  obtain,  as  shown  in  ftirt  II, 


E[Z^(t)]  ~  n^  Ae^°* 


where 


A  = 


i"(i)  r 

■>0 


dG(u) 


1-m  /  e"’^*^dG(u) 


/ 

Jo 


Define 


Then 


Gj^(si,S2,t)  =  2^  P[Z^(t)=Ji  Z^(t)=k]s^S2  • 

j^k=0 


Gj^(Bi,S2,t}  =  s^P^(t)  +  S2Pg(t)+l-G(t)-P3^(t)-Pj^(t) 
+  ^  li(Gp(s^,B2,t-u))dG(u)  . 


Since 


^  G”— (s-  j  S-^  *fc) 


*1  ""“s 


Si=S2=l 


we  obtain 


E[Z^(t)Zp(t)]  ~  n^rij^e 
This  suffices  to  show  that 


2cxb 


lim  p(W^(t),W  (t))  ^  1 

t->oo 
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Define 


G  (s^^Sg^tyr)  =  J  P[2;j^(t)=J;  .Zp(t+T)=k]s^s| 


d,k=-G 


Then 


Gj^Cs^^Sg^t^r)  =  S^SgPpCt+t)  +  l-CJ(t+T)-Pj^(t+T) 

+J  h(Gj^(sj^^S2,t-u,T))dG(u) 

rt+T 

J  h(Gj^(l,B2,t+T-u^0))dG(u)  . 


+  s. 


Since 


2 

^  Sg;t^  t) 


Si=8gfl 


we  obtain 


E[Zj^(t)Zp(t+T)]  -  Ae^^‘*'“'^ 


which  suffices  to  show  that 


lim  E[W„(t)-W„(t+T)]^  =  0 

t-^M  “  ^ 

and  hence  ^("t)  ^  ■»  variable,  in  quadratic  mean.  Also 

Var[Wj^]  =  A-1  >  0,  by  [3],  Ch.  6,  sect.  19,  so  that  is  non¬ 

degenerate  . 

That 

Ito^  p(W^(t),^^{t))  =  p(W^,Wj^)  =  1 

may  be  obtained  by  the  methods  of  Part  II.  Since  E[¥„]  =  1,  E[V??]  =  A, 

K  K 

a.e.,  R=l,2, ..,,n. 
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PART  II  -  TOTAL  NUMBER  OF  BIRTHS,  N(t) 


The  basic  model  may  be  any  of  those  previously  considered,  since 
we  will  now  work  only  with  the  lifetime  distribution  function  of  the 
entire  cell,  denoted  as  before  by  G(t)  =  P[cell  lives  for  age  <  t]. 
We  wish  to  study  the  total  number  of  cells  born  by  t,  denoted  by 
W(t).  See  Kendall  [4],  [5]  for  early  results  in  special  cases. 

To  facilitate  computations,  we  Introduce  the  following  variation 
on  the  basic  model,  as  indicated  by  the  figure. 


A  proliferating  cell  of  the  type  considered  in  previous  models, 
denoted  "type  2",  with  life  distribution  G(t),  divides  into  other 
"type  2"  cells,  whose  niamber  is  determined  by  the  generating  function 
h(s),  h'(l)  =  m  >  1.  However,  instead  of  considering  that  the  original 
cell  disappears  on  division  as  before,  we  suppose  it  is  replaced  by  one 
"type  1"  cell  of  infinite  lifetime.  Throu^out  we  assume  that  at 
t  =  0  we  initiate  the  cell  growth  process  with  one  "type  2"  cell. 

All  generating  functions  set  up  are  conditioned  on  this  fact. 

2.1  m(t),  Var[E(t)] 

Let  Z^(t)  =  number  of  type  1  cells  at  t 

ZgCt)  =  number  of  type  2  cells  at  t  . 
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Then  Z^(t)+Zg(t)  =  N(t)  =  total  number  of  type  2  cell  births  by  t. 


in^(t)  =  E[Z^(t)]  ,  mgU)  =  ECZgCt)]  ,  m(t)  =  E[N(t)]  . 


Define 


Then 


uo 

GCs^^Sg^t)  =  2  (t)=J}  Zp(t)=k]s^Sp  . 

J,k=0 


ft 

GCsi^Sg^t)  =  SgCl-GCt))  ■*■  SjL  J  ii(G(s^,Sg,t-u) )dG(u)  . 


Since 


SaCs^jSg^t) 

- ST - 


Si=S2=l 


=  m^(t) 


and 


3G(s^,Sgj t) 


Bi=S2=1 


=  »,(*)> 


adding^ 


iii^(t)  =  G(t)  +  mj  m^(t-u)dG(u) 

ra^Ct)  =  l-G(t)  +  mJ  m2(t-u)dG(u) 

m(t)  =  1  +  m  f  m(t-u)dG(u)  . 


■/ 

•^0 

Assuming  h’  (l)  s  m  >  l,  by  lermm  1  of  the  appendix,  for  t  large 

f  G(u)e"^V 

Jo 


m^(t) 


m  f  ue~^%G(u) 

_  Jo 


Qt 


m^a  f  ue”*^^  dG(u) 

Jo 


=  n^e 


aig(t)  ~ 


f  [l-G(u)  ]e"“% 

Jo 


mj 

Jo 


ue"^%G(u) 


Ob 


at 


/  -czu 
m  a  /  ue 


Jo 


dG(u) 


=  n^e 


at 


m(t)  ~ 


T 

Jo 


-au, 

e  du 


m  f  ue"^^dG(u) 

Jo 


at 


.at 


ma  /  ue”^dG(u) 


CO 

a/  I 

Jo 


=  n^e 


at 


and 


”l  "2  '  ”0 


For  h"(l)  <  w.  we  know  ([5]^  Ch.  6,  sect.  19^  21)  that 


2p(t) 

n^e 


a  random  variable,  in  quadratic  mean.  In  addition,  if 


f 

Jo 


E[W2(t)-Wg]  dt  <00  , 


then  WgCt)  almost  surely.  We  prove  a  corresponding  result  for 

W(t). 


2.2  Convergence  of  W.(t). 

”  0 

Theorem.  If  h'  (l)  =  m  >  1,  h"(l)  <  oo 


„  (t) .  , 

'  at  ”o  ' 


n  e 
o 
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a  random  variaTale,  in  quadratic  mean* 

Rroof .  Define 

CO 

F(s,,s„,t,'r)  =  J  P[N(t)=j;  K(t+T)=k]sW 
^  0<J<k 


Then 


rt  r t+T 

F(s^,S2,t,T)  =  B^SgCj  h(F(B^,B2,t-u,T))dG(u)+J  hCKCsg^t+T-u)) 

dG(u)+l-G(t+T  )] 


where 


00 

K(s,t)  ^  J  P[N(t)=j]s'^  . 

J=1 


¥e  wish  to  find 


3"^F(s^,S2,t,T) 


M2(t,T)  =  E[N(t)N(t+T)]  =*”5^-^ 


1  “^2 


=1==2=1 


Performing  the  indicated  differentiations, 

t+T  /•  t 


r  t+T  /- 1 

I^(t,T)  =  1  +  mj  m(t+T--u)dG(u)  +  mj  m(t-u)dG(u) 

T)dG(u)  . 

Using  the  method  of  proof  of  theorem  l8.1  of  Ch.  6,  for  t  large. 


+ 


h"(l)J'  m(t+T-u)m(t-u)dG(u)  +  nj'  I^(t-u,‘ 


M2(t,T) 


n^h" (1) 
o  ' 

lo 

1 

H 

1 

e"^°^dG{u) 

D 

2cit+cirr 


Now  we  can  compute 
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lira  E[W^(t+T)-W^(t)]  =  lira 


CO  ^ 


i_  pfl^Ct+T)  _  2Nft)H(t+T)  N^(t) 
^2  [^2a(t+T)  "  "^a(t+T)+at  2at  , 


h"(i) f 
Jo 


dG(u) 


l~mf  e”^^^dG(u) 

Jo 


1  +  1-2 


^2at+aT 

2cri;+0'T 


=  0,  as 


r  e"^“^G(u)  <  f  e"“%G(u)  =  i  . 
Jo  Jo  ® 


Hence j  completeness  shows  that 


W^(t)  =HitL 

O  Qt  O 

n  e 
o 


In  quadratic  mean.  The  theorem  is  proved. 

The  first  two  moments  of  W  are  obtained. 

o 


E[W  ]  =  lim  E 
^  t-»  00 


t-»eo  n  e  J  , 

o  1 


fshLl  - 1 

ne“* 

o 

h"  (1)  f  e^^^'^GCu) 
^  -^0 _ 

1-m  f  e  ^^^G(u) 

Jo 


VartW  ]  = 
o 


[h''(l)+m]r 

_ Jo _ 

1-m  f 

Jo 


dG(u)-l 


0  by  theorem  19-1  of  [3]^ 
Ch.  6. 


For  completeness  we  state  theorem  21.1  of  [5],  Ch.  6,  and  its 
corollary,  which  hold  here. 
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Theorem.  For  m  >  1,  h''(l)  <  m 


> 


G  not  a  lattice  distrihution,  and 


/: 


E[w^(t)-W^]  dt  <  00  , 


then  almost  surely. 

Corollary.  For  m  >  1,  h"(l)  <  o>,  G'  (t)  =  g(t)  and  such  that 


roo 

(  (g(t))^dt  <00  for  some  p  >  1  , 

Jo 

then  W  (t)  -^W  almost  surely. 


2,5  Correlation  between  W  and 

■■  ■"  '  "  o  c 

¥e  now  wish  to  find  the  correlation  between  the  random  variables 

W  and  W_ . 
o  2 

Recall  that 


G(sT,Bp,t)  =  2  PLz.(t)=^}  Zg(t)=k]s^s| 

^  J,k=0 


satisfies 


G(B^,S2,t)=S2(l-G(t))  +  s^^  h(G{s^,S2,t-u)dG(u))  . 


Since 


^  GCs^^Sg/t) 


H2i(t)  =  »  E[Z^(t)Zg(t)]  =  — 


1  “"2 


Bi=S2=1 


p  S‘^G(s  ,s„,t) 
M^^(t)  =  E[Z^(t)]^  - - - 


5s: 


Si=S2=l 


5  G(sT,s„,t) 


M22(t)  =  E[Zg(t)]^-  =  2 

5s„ 


s^=S2=l 
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we  obtain  the  following  integral  eguations. 


^  r't  p  rt 

^l^t)  =  2mj^  in^(t-u)dG(u)  +  h"(l) Tn^(t-u)dG(u)  +  mj  t^^(t-u)dG(u) 

in2(t-u)dG(u)  +  mJ'  I^2('t-u)dG(u) 

r  mj^(t-u)in2(t-u)dG(u) 


•f 


u)dG(u)  . 


The  method  of  theorem  l8.1  of  [5],  Ch.  6,  yields,  for  t  large. 


h"(l)  r^Te-^iaM 
_ ^ _ 

1-mf  e"^^dG(u) 

Jo 


h"(l)  re"^“'^G(u) 

I  1-m  /  e  ^^dG(u) 

Jo 


Mon (t) 


h'’(l)n^n^J  e'^'^dGCu) 

l-mre-2«V(u) 

Jo 


Since  W(t)  =  Z^(t)  +  ZgCt)  , 


59 


p(t)  H  p(Z2{t),N(t))  = 


E[Z^(t)+Zg(t)  ]  [Z,,(t)  ]-E[Z^(t)+Z2(t)  jELZgCt)  ] 

] 


-*2^ «/ 

/(E[W^(t) ]-E^[N(t) ]) (E[Z^(t) ]-E^[Z.(t) : 


Let 


then 


For 


K(s,t)  -  2  P[W(t)=j]s'5 
J=1 


K(s,t)  =  s[l-G(t)  +  h(K(s,t-u))dG(u)]  . 

Jo 


l^(t)  =  E[K^(t)]  = 


s=l 


,  we  obtain 


ft^(t)  =  2mJ^  m(t-u)dG(u)  +  h"(l) J  m^(t-u)dG(u) 
+  m^  l^(t-u}dG(u)  . 


For  t  large 


t^(t)  - 


n^"  (1)  f  e"^“%G(u) 

°  Jo 


Jo 


1-m  j  e^^HiG(u) 


2at 


so  that  writing 


p(t) 


M^g(t)  +  -  m^(t)ni2(t)  -  v^it) 


Putting  c  =  /  e“^^%G(u) 


-JO 


h"(l)n^n2C  +  h" (l)n2C 


1-mc 


2 

n^^ng^n^ 


iA"(i)<=  g 


1-mc  o 


'^n|h"(l)c  2' 

1-mc  ”  ^2 


itO 


11m  p(t) 

X— >  CO 


since  we  have  assumed  all  along  that  G  is  not  a  lattice  distribution, 

we  may  factor  out  -l)  >  0  in  the  expression  for  p(t).  See 

Ch.  6,  sect.  19. 

Thus 


pCt) 


ng(n^+ng)(  -l) 


=  1 


We  wish  to  conclude  that 


Coviy  (t),W  (t)]  Cov[W.¥„] 

p(t)  - - 2 - _  - 2-^  =  1  . 

cr[W^(t)]orW2(t)]  c[W^]a[W2] 

Since 

E[W^(t)-W^]^  ^0  , 

E[W2(t)-W2]^  -»0  , 

we  know  that 

E[W^(t)]  -^E[W^']  , 

E[W2(t)]  ^£[¥2]  . 

Then  it  suffices  to  show  that 

|E[W^(t)W2(t)]  -  E[W^W2]|  ->0  . 

This  follows  from 

Wo(t)W2(t)-W^W2  =  W^(t)W2(t)-W^W2(t)-tW^W2(t)-W^W2 

=  W2(t)(w^(t)-WQ)  +  WQ(W2(t)-W2)  . 


We  have  by  the  Cauchy-Schwarz  inequality 


IeC’W  (t)¥2(t)--W^W2]|^  <  E[V^(t)]E[W^(t)-W^]^+E[w^]E[W2(t)-W2]^  -^0 
since 

E[W^]  <w  and  E[Wg(t)]  -»E[W2]  <  “  . 

Hence 

PCW^.Wg)  =  1  . 

Comparison  of  the  first  two  previously  derived  moments  of 
with  the  corresponding  moments  of  Wg  given  in  [3]^  Ch.  6,  sect.  19 
shows  that  they  are  the  same. 

Hence 

Wg  =  almost  surely. 

2.k  m(n)  in  discrete  case. 

To  check  some  of  the  results  for  the  means  m^Ct)  and  mgCt), 
suppose  we  consider  the  corresponding  discrete  case,  in  which  a  cell 
divides  at  the  end  of  a  unit  interval  of  time. 

f  V12I 

The  expected  value  matrix  M  =  \^2Jl^22\  ’ 

m^  ^  =  E[number  of  cells  of  type  j  issuing  from  a  type  i  cell 

at  its  mitosis] 

m  /  . 


for  ^,1=1,2 


M  = 


fl  0  1 

1  m 


I  1 
m  -1 


m-1 


k2 


> 


and  define  a  >  0; 


Let  G(u)  =  U(u-l)^  where  U(r)  = 


'1,  r  >  0 
^0,  r  <  0 


X' 


e’^%G(u)  =  e 


-a  ^  1 
”  m 


and  we  obtain 


n  on 
m  =  e 


Substituting  G(u)  =  U(u-l)  into  expressions  for  m^(t),  m2(t) 
previously  derived,  and  changing  t  to  n,  n  large. 


an 


”*1^*^^  ~  m  ln(m) 


/  \  (m-l)  an 

“a'"'  ~  ® 


from  continuous  case 


nigCn) 

m^Iriy 


~  m-l,  which  checks  with  the  matrix  case,  in  -vdiich 


an 


”i<"’  ~  S3 


m^(n)  ~ 


an 


2.5  m(t)  for  m=l,  m  <  1. 
Let 


K(s,t)  =  t  PtN(t)-d]s'^ 
J=0 


K{s,t)  =  s[{l 


-G(t)  +[ 
-^0 


h(F(B,t-u) )dG(u) ] 


m 


(t)  -  E[N(t)]  =  Msiii 

OS 


s=l 


For  the  case  m=l, 

ia(t)  =  1  +  f  itt(t-u)dGCu) 

•Jo 

and  the  system  Is  mathematically  equivalent  to  a  renewal  process.  For 
mild  restrictions  on  G,  such  as  /J*  tdG(t)  <  ^  G(0)  =  0, 


m(t)  ~  . -  ■"  ,  t  large  . 

[  udG(u) 

Jo 

See  Chapter  ?  of  Bellman  and  Cooke  [1]. 

Case  m  <  1. 

m(t)  =  1  +  m  (  m(t-u)dG(u)  . 

Jo 


By  lemma  2  of  the  appendix,  for  t^«> 


m(t)  -> 


1 

1-m 


We  check  these  results  for  cases  m=l,  m  <  1,  by  reverting  to  the 
discrete  case. 


2.6  Convergence  and  moments  of  N(t),  m  <  1. 

We  briefly  treat  the  case  m  <  1. 

Theorem.  Let  G  be  a  distribution  function, 

G(0)  =  0. 

For  m  <  1, 

N(t)f  W^,  a  random  variable,  a.e. 


1^4 


Proof.  NCt)^  a.e. 


Since  E[N(t)J  =  m(t)  satisfies 

ni(t)  =  1  +  m  /  m(t-u)dG(u) 

JO 

m(t)  ss  t  -♦<»,  "by  lemma  2  of  the  appendix. 

Hence  N(t)  is  finite  a.e.  for  all  t.  After  that  point  t 
which  ZgCt)  =  0  a.e., 

N(t)  =  N(t+T)  <  00  for  all  T  >  0  a.e. 


Hence  N(t)'f  a  random  variable  a.e. 


E[1J  ]  =  ^  . 
o  1-m 


To  obtain  hi^er  moments  of  M  we  assume  h^“^(l)  <  oo  for 

o 

2^  •  •  •  » 


Define 


K(s,t)  =  f  P[W(t)-j]s'^  . 
J=1 


From 


K(s,t)  =  s[l-G(t)  +1  h(K(s,t-u))dG(u)l 

Jo 

we  obtain  for  W^(t)  =  E[K^(t)]  the  equation 

l^(t)  =  2m(t)-l  +  h"(l)J’  m^(t-u)dG(u)  ■*"  l^(t-u)dG(u) 


Claim: 


f  m^(t-u)dG(u)  -»  (-5-^)^  as  t->«> 
Jo  -L-m 


for 
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Proof  of  claim:  For  0  <  7  <1, 


^  m^Ct-u)dGCu)  +/  ^  m^(t-u)dG(u) 

Jt’' 


r  m^(t-u)dG(u)  =f 
Jo  JO 

Since  m(t) 'f  as  ^ 


J'^m^(t-u)dG(u)  <  ra^(t)  [GCt)-G(t^)]  ^  (^j^)^[G(t)-G(t’')  ]  -->0  , 


and 


m 


r^7 

■(t-t^)G(t^)  <1  m^(t-u)dG(u)  <  m^(t)G(t^) 

"Jo  ” 


\diich  suffices  to  establish  the  claim. 
Hence 


where 


aM 


j^(t)  =  A(t)  + J  f^(t-u)aG(u) 

BB  t-. 


^  1-m  [l-m 


o 


Hi^er  moments  are  similarly  obtained  recursively,  but  no  general 
form  for  the  generating  function  of  has  been  obtained. 


2.7  m(n)  for  m=l,  m  <  1,  in  discrete  case. 
For  m=l 

fl  0\  .n  /I  0\ 


and 


m(n)  =  n+1  n  for  n  large. 


For  m  <  1 


M  = 


(1  0. 


0 


n 


m 


/ 


m(n)  = 


n 


m 


+  m' 


n 


1 

l-m  ^ 


as  ri->  w 


which  checks  with  the  continuous  case. 


2.8  Moments  of  N(t),  case  m=l. 

We  consider  the  asymptotic  moments  of  N(t)  for  a  simple  age 
dependent  branching  process  with  m=l.  Two  theorems  from  Laplace 
transform  theory  will  be  useful.  See  Widder  [10]. 

Abelian  theorem.  If  for  some  y  >  0 


and 


lla  eM  =_£_ 
t-»M  ^7  r(y+l) 


then 


*•« 

li(s)  =  e”^^dm(t) 

Jo 


11m  s^pi(s)  =  c 
siO 


Tauberian  theorem.  If  mCt)^  is  such  that 


r”  -St 

ia(s)  =  I  e"  dm{t)  converges  for  iR^Cs)  >  0 

Jo 
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and  If  for  some  y  >  0 


lim  s^|j,(b)  =  c 
a^O 


then 


t-»  0#  r  (7+1) 


Lemma.  Let  G  he  a  distrihutlon  function,  G(0)  =  0 
Let  H(t)t  satisfy 


H(t)  =  f(t)  +f  H(t-u)dG(u) 

Jo 


■vdiere 


f  (t)  is  bounded  on  every  finite  interval,  and 


lim 


&  • 


Then 


lim 

w 


H(t) 


n+l 


T^hoT 


ra, 


G 


Proof.  By  [5],  Gh.  6,  appendix,  the  unique  solution  bounded  on  every 
finite  interval  is 


■where 


f(t-u)dM(u) 


M(u)  =  2  G^“^(u)  . 
n=0 


Denoting  the  Laplace  transform  of  a  function  g  by  g,  ■we  obtain 
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1 


H(s)  =  ?(s)  ll(s)  =  f(B)  ^ 

Ii-g(b)  J 

As  s^O  ,  a  Taylor  expansion  to  the  first  moment  yields 

I-GCs)  -  sm^ 

and  hy  the  Ahelian  theorem,  as  s^O, 


?(b)  ~  sSslil 

‘  n 


H(s)  ~ 


ar(n+l) 
n+1 


m, 


‘G 


and  hy  the  Tauherian  theorem 


11m  H(t)  ar (n+l)  _  a 

t-^co  ^n+1  (n+2)  “  (n+l)mQ^ 


We  may  now  obtain  asymptotic  moments  of  W(t)  for  m=l. 
Theorem.  If  m=l. 


h(i^)(l)  <e,  ^  n=l,2,  ... 


/; 


udG(u)  H  m.Q  <  00  , 


k9 


then 


£[11^(1})]  —  M^Ct)  >  ••• 


M^(t)  a^[h''(l)]‘ 


•t-»»  •(; 


2n-l 


2n-l 


«here  aj^=l,  ag=  3,  b.^  3,  a^  =  ^,  a^  =  ^,  and  the  (a^)  may  be 


obtained  recursively. 


Proof.  Define 


K(s,t)  =  X  P[N(t)=jls'^ 

J=1 

rt 

K(B,t)  =  s[l-G(t)  +  I  h(K(s,t-u))dG(u)] 

•^0 

M^(t)  =  m(t)  =  1  +  J  m(t-u)dG(u)  . 


From  Bellman  and  Cooke  [l],  pp.  236-239/ 


^  “g 


MgCt)  =  2[m(t)-l]  +  h"(l) JV(t-u)dG(u)  +J^  M2(t-u)dG(u) 

f\^(t-u)dG(u)  =  f  m^(t-u)dG(u)  +f  m^(t-u)dG(u)  for  0  <  y  <  1, 

Jo  Jo  Jt^ 

N(t)‘f  ,  hence  and  ,  and 

f  m^(t-u)dG(u)  <  m^(t-t^)  [G(t)-G(t’')]  =  o(t^) 


m^(t-t^)G(t^) 


2 


m  (t-u)dG(u)  <m  (t)G(V) 
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BO  that 


Mgft)  =  A(t) 


idiere 


By  the  lemma. 


lim 

t-^  03 


A(t)  h'Kl) 

.2  "  2 

t 


lim 
t— ^  00 


^2^^^  =  h” (1) 
t^  5 

G 


The  general  result  follows  recursively  in  the  same  way.  No 

general  formula  for  the  {a^]  has  "been  obtained. 

Note  that  the  result  for  I^(t)  differs  markedly  from  the  renewal 

theory  case,  where  h(s)  =  s,  and  with  the  second  moment  asymptotic  to 
2 

t  Eind  the  variance  asymptotic  to  t. 


2.9  Moments  of  N(n)  in  discrete  case,  for  m=l,  m  <  1. 

The  two  theorems  to  he  given  in  this  section  are  limit  theorems 
concerning  the  total  number  of  births  in  discrete  time  for  cases 
m=l  and  m  <  1.  These  results  are  similar  to  those  just  given  for 
continuous  time. 

A  cell  living  at  time  n,  n=l, 2, . , .  will  devide  at  time  n+0 

into  K  cells  with  probability  p„,  K=0,l,2, ...  and  we  let 

A 

00  /  \ 
h(s)  =  ^  f  h'(l)  =  m,  and  suppose  h^'^^(l)  <  » 

k=C  ^ 

for  i=l,2, ...  and  h"(l)  >  0. 
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Independence  as  before  Is  assumed.  Each  new  cell  will  divide  In 
accord  with  h  at  time  (n+l)  +  0. 

Let  lT(n)  be  the  total  number  of  births  by  time  n,  n=l,2, ... 
starting  with  one  cell  at  time  0. 

Define 

00 

K(s,n)  -  2  P[N(n)=j]s'5  ,  n=l,2,... 

J=1 


By  taking  G(u) 
time  t,  t  >  1 


1,  u  >  1 
p,  u  <  1 


in  the  equation  in  continuous 


h(K(B,t-u))dG(u)] 


and  setting  t=n+l,  we  obtain 


K(s,n+l)  =  s  h(K(s,n))  . 


We  proceed  to  cases  m=l  and  m  <  1. 

To  treat  the  case  m=l,  we  need  the  following  particular  case  of 
the  lemma  of  Ifert  II  applied  to  discrete  time. 

Corollary  1.  Let  ff  be  a  function  on  the  positive  integers 
satisfying,  for  K  n=l,2, and  a  >  0  and 

f(K+l)  =  g(K)  +  f(K) 

lim  =  a  . 

K-^co 


Then,  as  K  ->«, 


Theorem  1.  Let 


E[N(n)  ]  =  M(n) 

E[lf(n)]  =  M^(n)  ,  r=2,5,... 

Then  If  m=l,  we  have 

M{n)  =  n 

and 

M  (n)  ^  , 

lim  =  a^[h"(l)] 

where  the  [a^3^_2  ^re  the  a'e  in  the  analogous  continuous  time 
theorem  for  If(t),  m=l. 

Proof  sketch.  From  K(s,n+l)  =  Bh(K(s,n)),  since 

^K(s,n)  I  =  M(n)  , 

B=1 


Hence 


M(n+1)  =  1  +  M(n)  ,  m{l)  =  1 

M(n)  =  n 


M2(n+l)  =  2  M(n)  +  h"(l)  M^(n)  +1^2(11)  • 

By  the  lemma,  since  N(n)t  implies  I]^(n)|  implies  I^(n)T  r 


lim 

n-j  00 


k — 1,2, #  # •  « 


Continuing  in  this  recursive  manner,  the  result  may  he  estahlished. 

The  results  for  the  case  m  <  1  are  essentially  identical  with 
those  for  the  continuous  case. 
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Theorem  2. 


If 


m  <  1,  as  n 


M(n) 


1 

1-m 


M^(n)  r=2,3, ... 

00 

where  the  <jh  >  «  are  the  same  as  in  the  continuous  case,  m  <  1. 
r  r=c 

Proof.  From  K(B,n+l)  =  sh(K(s,n))  , 


M(n+l)  =  1  +  m  M(n)  . 

We  may  apply  lemma  2  of  the  appendix  for  G(u) 
to  obtain 


■{ 


1 ,  u  >  1 
0,  u  <  1 


M(n)  . 


J^(n+l)  =  2mM(n)  +  h"(l)M^(n)  +  mJ^(n) 


and  again  applying  lemma  2  of  the  appendix. 


-m 


and  so  on  reciirsively  to  establish  the  result. 


1+m  h'*(l) 

(l-m)® 


S'* 


BLET  III  -  TWO  TYEEJS  OF  CELLS 


3.1  Model  I. 

Tills  section  discusses  two  models  for  cell  growth  In  an  expanding 
population  involving  the  proliferation  of  two  types  of  cells.  In  the 
first  modelj  a  "type  1"  cell  has  a  non- lattice  life  distribution  H, 
H(0)  =  0  and  is  incapable  of  proliferation.  At  the  end  of  its  life 
period,  it  is  absorbed.  A  "type  2"  cell,  with  non-lattice  distribu¬ 
tion  G,  G(0)  =  0  is  capable  of  proliferation.  At  the  end  of  mitosis 
it  gives  rise  to 

2 

2  clusters  of  "type  1"  cells  -  with  probability  q 

1  cluster  of  "type  1"  cells  and  a  cluster  of  "type  2"  cells  -  with 

probability  2pq 

2  clusters  of  "type  2"  cells  -  with  probability  p^ 

where  p+q  =1,  p  >  0,  q  >  0,  and  a  cluster  of  type  2  cells  is  a 
collection  of  those  cells,  containing  k  cells  with  probability  p  , 
k=0,l, ...  .  Define 

00  , 

h(s)  =  2  Pv®  f  h' (l)  =  m,  h"(l)  <00  . 

k=0  ^ 

Similarly,  a  cluster  of  type  1  cells  contains  k  type  1  cells  with 
probability  q^^,  and  define 

■b(s)  =  2  q,6^  ,  b*  (1)  =  d 
k=0  ^ 

Each  type  of  cell  and  cluster  forms  independently  of  the  state  of  the 
system.  The  cells  are  independent  of  each  other.  See  Kendall  [6]. 
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For  h(s)  =  bCs)  =  s  the  process  may  he  represented  schematically 
as  follows 

2  o  2 

P  2pg.  q 

Define,  given  that  at  t=0,  the  process  starts  with  one  newly 
horn  type  2  cell, 

Z^(t)  =  number  of  "type  1"  cells  at  t 

ZgCt)  =  number  of  "type  2"  cells  at  t 

NgCt)  =  nuniber  of  "type  2"  cells  horn  hy  t 

N^(t)  =  niuribei’  of  "type  1"  cells  horn  hy  t 

5.2  E[Z^(t)],  ELZgCt)]  . 

The  first  moments  of  these  quantities  are  obtained  for  large  t, 
and  are  related  to  the  corresponding  discrete  time  case  results.  A 
discussion  of  convergence  in  quadratic  mean  Indicates  also  the  asymp¬ 
totic  second  moments. 

Before  starting  the  computations,  we  note  that  the  "type  2"  cells 

form  a  process  which  is  independent  of  the  type  1  cells  as  follows. 

2 

With  prohahility  q  -  0  "type  2"  cells  are  emitted  on 

division  of  a  "type  2"  cell 

2pq  -  1  cluster  of  "type  2"  cells  is  emitted 
2 

p  -  2  clusters  of  "type  2"  cells  are  emitted. 

The  equivalent  generating  function  for  the  numhei*  of  new  "type  2" 
cells  created  on  division  of  a  "type  2"  cell  is 
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f(r)  =  g(h(r)) 


2  2  2 

where  g(s)  =  q  +  Spqs  +  p  s  and  h(r)  =  2  • 

J=0  -J 

Hence,  as  done  in  previous  sections,  the  backward,  forward,  and 
total  time  distributions  and  asymptotic  mean  and  variance  of  the  number 
of  "type  2"  cells  in  a  state  may  be  computed,  if  the  type  2  cell  pro¬ 
gresses  throu^  n  states,  the  state  life  characterised  by  its  dis¬ 
tribution  F,,  k=l,2,  ...,n.  Here  G  =  F.  *»••*  F  if  we  deal  with  the 

XU 

sequence  of  states  model,  and  G  is  the  equivalent  life  distribution 
of  the  cell  in  the  semi -Markov  case. 

Throughout,  all  generating  functions  depending  on  the  time  are  con¬ 
ditioned  on  the  event  that  at  t=0,  the  process  begins  with  one  newly 
born  cell  of  type  2  or  one  of  type  1,  as  will  be  specified. 

Define 

G  (8  ,8  ,t)  =  f  P[Z,  (t)=J)  Z„(t)=k|''yS!_J  “*]8W 
^  ^  J,k=0  t-u  X  ^ 


Let 


Gj^(6i,S2,t)  =  s^[l-H(t)]  +  H(t) 

G2(Si,S2,t)  =  B2(l-G(t))  +J^  [q^^(G^(B^,S2,t-u))+2pqb(G^(s^,S2,t-u)) 

h(G2(Bi,Bg,t-u))  +  p^^(G2(s^,S2jt-u))]dG(u) 

^G2(B^,B2>t) 


m^(t)  s  E[Z^(t)]  - - ^ 


s^=S2=l 


SGo(b  ,Sp,t) 

m2(t)  =  ELZgCt)]  =  — - - 


Si=E2-1 
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or 


in^(t)  =  2qdJ  [1- H(t-u)  ]dG{u)  +  2pm J"  in^(t-u)dG(u) 
Tn^(t)  =  2qd[G(t)  -  G  *  H(t)]  +  2paiJ  iiu^Ct-u)dG(u) 


m^Ct)  =  l-G(t)  +  2pm  J  m^(t-u)dG(u)  . 


f 

^0 


To  treat  the  case  of  an  expanding  cell  population,  assume  that  2pm  >  1, 
and  we  define 


r 

a  : 

Jo 


e-'^dG(u)  = 


2pm  • 


Then,  for  t  large,  hy  lemma  1  of  the  appendix, 


m^(t) 


I 


2qd  /  [G(u)-(G*H)  (u)  ]e  ^du 


f 

Jo 


2pm  /  ue""^'^  dG(u) 


at  at 
e  =n^e 


r 

Jo 


[l-G(u)  ]e"°^^du 


2pm  /  ue”^^  dG(u) 


i 

Jo 


at 


i[i_ 

2pm^^ 


2pm 


f  ue  ^%G(u) 

Jo 


=  n^e 


at 


3.3  Comparison  with  discrete  case. 

We  compare  these  results  with  those  in  the  discrete  case- 
Let  X  =  number  of  "type  2"  cells  emitted  from  a  "type  2"  cell  on 
division. 

E[X]  =  E[xj^  type^l^ ^cluster j  pj-^  "-type  1"  cluster  emitted]  + 

E[x|0  "type  1"  clusters  emitted] 

P[0  "type  1"  clusters  emitted] 

E[X]  =  2pqm  +  p^2m  =  2pm. 
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Let  Y  =  number  of  "type  1"  cells  emitted  \dien  a  "type  2"  divides. 

E[Y]  =  d.2p(i  +  2d.q^  =  2qd  . 

The  expectation  matrix  M  is 

^  “  2qd  2pm^ 

■vdiere  we  assume  for  simplicity  that  H  =  0,  that  is,  a  tjrpe  1  cell  has 
Infinite  life. 


= 


Let  G(t)  =  U(t-l)  ,  U(s)  = 


1,  B  >  0 


_0,  B  <  0 

Let  t  correspond  to  n,  the  n^^  step  of  a  discrete  process.  Thus 
the  continuous  case  is  changed  to  the  discrete . 


-a 


2pm  ' 


m^(t)  ~  m^(n 


2qd 

)  -  e™  =  .3^  e““ 

'  1  pna 


m^Ct)  ~  i^(n)  ~  i  [1-  g^le“" 


"l*")  ^  2ad_ 

nigtn)  ~  2pm-l 


From  the  matrix  M, 


~  2i^  ^ 


an 


2pm' 

m^Cn)  ~  (2pm)’^  =  e^“ 
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The  ratios  check , 


m^TnJ  2pm-l 

This  leads  to 


3.^  Comparison  vlth  results  of  Snov. 

Iiemma.  The  ratio  of  the  components  of  the  left  eigenvector  cor¬ 
responding  to  the  largest  eigenvalue  of  M  is  that  of  the  asymptotic 
ratio  of  the  mean  numbers  of  "type  1"  and  "type  2"  cells  in  the  contin¬ 
uous  case,  in  T^ilch  H  =  0  [i.e,,  type  1  cells  have  infinite  life]. 

That  is 


if 

then 

Eroof . 


a 

p  ~  ^  ^  • 


1  0 

^  "  ^2q,d  2pm^ 


X 

max 


=  2pm 


(a,p)M  =  2pm(a,P) 


q:  +  2g,dp  -  Cjsao: 

2  - 

p  2pm-l  ~  mg(t) 

Snow  [8]  presumably  has  general  results  of  the  type  in  the  lemma, 
but  for  the  "irreducible"  case,  that  is,  if  M  is  now  normalized  to 
form  a  transition  matrix,  it  would  be  irreducible.  In  the  case  treated 
here,  the  states  of  M  do  not  communicate.  See  also  [3]j  Ch.  5* 
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To  find  E[N^(t)],  set  E  =  0,  so  that  Sg^t)  =  -which 

is  equivalent  to  letting  a  type  1  cell  have  infinite  life. 

E[N^(t)]  ~ 

Ihese  results  can  he  checked  vith  the  discrete  case^  and  a  lemma 
like  the  above  proved. 


6l 


3.6  Convergence  in  quadratic  mean  and  asymptotic  correlations. 
Given  h"  (l)  <  <»,  2pm  >  we  wish  to  show  that,  t-»a), 


Z  (t) 

W^{t)  - - ^  In  quadratic  mean 

nie^ 


W^Ct)  = 


Zo(t) 


at  2 


w,i(t)  = 


N,(t) 


at 


w. 


01 


where 


V<")  = 


N2(t) 

at 

^^02® 


w 


02 


“  ‘  Jo  ^  • 


The  computation  of  all  new  generating  functions  to  follow  assumes 
that  at  t  =  0,  the  process  starts  with  one  cell  of  type  2. 

We  shall  let  the  notations  G^Cs^^Sg^t)  and  G^Cs^^Sgjt)  stand 
for  the  generating  functions  which  they  stood  for  in  the  preceedlng 
section  on  two  cell  types. 

Define 


00 


B(s-,Bp,t,T)  =  X  P[Z  (t)=J>  Z  (t+T)=k]s^s5 
^  j,k=0 


Then 
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B(B^,82,t,T)  «  1-G(t+T)  +  j  [(l^^(G^(s^A^'t-u))+2pq>(G^(s^,l,t-u)) 


2,  2 


h(B(s^,S2,t-u^T))  +  ph  (B(Sj^,82jt-u,T))]dG(u) 


Ijt+T-u))  +  2pq;b(G^(s2A^t+T-u)  )h(G2(s2,l,t+T-u)  ) 


t+t 


+  p^h^(G„(s«,l^t+T-u)) ]dG(u)  . 


Since 


d  B(s^,S2^t,T) 


^=1  ^2 


=  E[Z^(t)Z^(t+T)]  , 


=2  »2 


we  find,  using  the  method  in  theorem  l8.1  of  [3]}  Ch.  6,  that,  for 
large  t, 

r 

2. 


Let 


E[Z^(t)Z^(t+T)] 


2pn^  f  e”^*^^G(u)  [h"  (l)+pm'^ 
■^-^0 


] 


r 


l-2pmj|^  e"^^%G(u) 


2at+aT 


C(6^,s  ,t,T)  =  2  (t)=ji  Zp(t+T)=k]s^S2 

^  J,k=0 

ft  2 

C(s^,B2,t,T)  =  s^S2[l-G(t+T)]  +j  [q  +2pq  h(C(s^,B2;.t-u,T)) 


2,  2 


+  p  h  (c(s^,S2>t-u,T))]dG(u) 


't+T 


2.„  ,  ,,  .  .  ...  2,2, 


+  s^j  [q  +2pqh(G2(l,S2,t+T-u))+p  h  (G2(l,S2t+T--u)  )  3dG(u) 


8  C(s,,Sp,t,T) 

- Sj^.S2=l  “  Etz^(t)Z2(t+T)  J 
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Sprig  f  e'^^^dGCu)  [h"(l)+pm^] 
E[Z2(t)Z2(t+T)]  ~  - 2 - - - gSorti+oT 

l-Spin  f  e**^^  dG(u) 

-^0 


QO 

D(s^,s„,t,T)  =  2  P[N^(t)=J;  NT(t+T)=k]s^Sp 

k>J>0-^ 

D(s^,S2,t,T)  =  1-G(t+T)  +J  +  Spqb(s^S2)h(D(s^,S2,t-u,T)) 


+  p  h  (D(s^,S2jt-u,T))]dG(u) 


J  ( S2)+2pql)  (s2)h(R(  Sg,  t+T-u) ) 


+  p  h  (RCsg, t+T-u) ) ]dG(u) 


vhere 


R(s,t)  =  D(s,l>t,D) 


a^D(s^, Sg^t^T) 


5s^  dSg 


=  E[N^(t)N^(t+T)] 


’1  “2 


E[N^Ct)N^(t+T)] 


SpriQ^  J"  e”^^%G(u)  [h”  (l)+pm‘^ 


l-Spn/  e”^^dG(u) 


Sat+cer 


6k 


00 


H(s^,S2 


where  we 


ft,7)  =  S^Sg 


f  t+T 

i 


way  take 


.t,T)  =  2  P[Np(t)=ji  Np(t+T)=k]B^sJ 

6  >1  ^  ^ 

ft  p 

[l-G(t+T)  +J  [q,  +2pqh(H(s^,  Sp^t-u^r) 

+  P^h^(H(s^,Sg,t-u,T))]dG(u) 
[q^+2pqh(K(B2,t+T-u)  )+p^h^K(s2,t+T-u) )  ]dG(u) 


K(s,t)  =  H(s,l,t,0) 


^H(Si, Sg^tjT) 


Ss^'Ssp 


=  E[N2(t)N2(t+T)] 

X  d 


E[N2(t)N2(t+T)  ~ 


2pn' 


02 


f 

Jo 


e  dG(u)  [h"  (l)+pin^] 


r" 

n 

Jo 


l-2pmf  e"^°^dG(u) 


2at+aT 


Applying  the  shove  formulas^  we  see  that 


E[W^(t)-W^(t+T)]^  ->  0 
E[W2(t)-W2(t+T)]^  ->0 
E[WQ3  (t)-WQ3^(t+T)]^  -^0 

EfWQ2(t)-WQ2(t+T)]^  ^0 

so  that  W^(t)  W^Ct)  -V  ‘l^dratic 

mean. 
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Ealrwlse  correlation  among  ,  Wg ,  . 

¥e  explicitly  indicate  how  to  find  the  correlation  only  het>reen 
¥^  and  ¥g.  The  other  five  follow  similarly. 


O  Gp(8^,Spjt) 


,  9  G„(s,,s„,t) 

M,  (t)  =  Et4(t)]  =  - - 

OBg 


Si=S2=l 


°  Sg^t) 


-  E[Z3^(t)Z2(t)]  =  - |s^=S2=l 


jJe-^^^dG(u)  =  c 

2pn^c  [h"  { 1 )  +pm^  ]  OfYt 

- J  - 

2pnpC  [h"  ( 1 )  +pm^  ]  2o± 

2pn  ngc [h" (1) +pm^  J  got 
^2^^^  ~  l-2pm.c  ^ 

M.p(t)  -  nL (t)mp{t) 

p(Z^(t),Z2(t))  =  ^  -■  -  — 

^(M^l(t)-ii^(t))  (M22(t)-m^{t)  ) 


p(Z^(t),Z2(t)) 


n^Hg  L  l-gpnic 


=  1  , 


n^ngl 


l-2pffic 
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since  ve  have  the  following  lemma. 


Lemma. 

2pc(h"(l)+pm^)  _  1  >  0  . 
l-2pmc 

Proof.  First  note,  hy  the  Cauchy-Schwarz  inequality 


=  f  e”^‘^%G(u)  >  [  [  e'’'^dG{u)  =  — i-? 
-^0  Jo  4pm 


2pc(h'’  (l)+pm^)  _  _  c  [2ph"  (l)+2p^m^+2pm]-l 

l-2pmc  ”  l-2pmc 


It  suffices  to  show 


c[2ph"(l)+2p^m^  +  2pm]  >  1 


or  to  shov 


or 


or 


Ph"(l) 


2  2 

+  2pm  + 


2  2 


2pm 


>  1 


h”(l)  +  pm  +  m 
2pra^ 


>  1 


h"(l)  +  m  >  pm^ 


hut 


O  p 

h"(l)  +  m>m  >pm  ,  as 
h'*(l)  is  the  second  factorial  moment  of  h. 


6? 


In  order  to  show  that 


Cov[W  (t),W-(t)]  CovCWj^W^] 

p(Wi(t),Wp(t))  - - - - - - » - =  1 

^  ^  a[w^(t)]a[W2(t)]  o[W3_}crrW2] 

we  use  the  same  reasoning  as  in  the  section  on  total  number  of  births 
with  one  cell  type. 

Since  E[Wj^(t)-W^]^  -^0 

->0 

we  have 

E[W^(t)]  ->E[W^] 

ElW^it)]  ->E[W2]  . 

It  suffices  to  show  that 

lE[¥^(t)W2(t)-W^W2]|  ->0  . 

Write 

W^(t)W2(t)-W^W2  =  W2(t)[Wj^(t)-Wj^]  +  W^[¥2(t)-¥2]  * 

By  the  Cauchy-Schwarz  Inequality, 

|E[W^(t)W2(t)]-E[¥^W2]|^  <  E[W2(t)]E[Wj^(t)-W;^]^  +  E[w^]E[w2(t)-W2]''-^  0 

since 

Et.^)  <  «  and  EEw^Ct)]  -^EfWg]  <  «•  • 

Hence 

=  1. 

All  the  other  five  correlations  are  1. 


% 
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since 


E[W^]  =  ECWg]  =  E[Wq^]  =  E[Wq2]  =  1 


and,  by  the  previous  section. 


Var[W^]  =  Var[W2]  =  Var[WQ^]  =  Var[WQ2]  = 


2p  r  e'^^^'^dOCu)  [h"  (l)+piii^+  m]-l 
■JQ _ 

l-2pin  f  e"^°^dG(u) 

Jo 


>  0 


Wi  =  Wg  =  =  Wo2 


a.s. 


In  the  discrete  case. 


m  >  1,  convergence  of  11m 


n->oo  X 


n 


-» w. 


k' 


max 


a  random  variable  in  quadratic  mean,  is  discussed  in  [3],  Ch.  2,  for 
multiple  types  k=l,2, . ,  .,k^.  Here  is  the  largest  eigenvalue 
of  a  relevant  expectation  matrix. 


5.7  Model  II. 

The  second  model  of  two  cell  types  to  be  considered  is  similar  to 
the  first.  Using  the  same  notation  as  in  the  previous  case,  the  pro¬ 
cess  may  be  described  as  follows.  At  the  end  of  mitosis  of  a  type  2 
cell,  the  progeny  may  be 

1  cluster  of  "type  1"  cells  and  1  cluster  of  "type  2" 

cells  -  with  probability  q 

2  clusters  of  "type  2"  cells  -  with  probability  p 

p+q=l. 
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In  case  h(s)=s,  l>(s)=s,  we  may  represent  the  process 


schematically  as 


5.8  E[Z^(t)],  ECZgCt)]. 


Let 


G  (s  ,B  ,t)  =  X  P[Z  (t)=J>  Z  (t)=k 


type  in  J  k 
at  t=0-'  1  2  ‘ 


Then 


G3^(si,B2,t)  =  s^[l-H(t)]  +  H(t) 

GgCs^^Sg^t)  =  B2[l-G(t)]  +J^  [ph^(G2(B^,S2,t-u)) 

+  qh(G^(s^,S2,t-u)  )h(G2(Sj^,B2,t-u))  ]dG(u) 
Assuming  an  expanding  cell  population,  or  m(p+l)  >  1, 

~  r“ 

qd  /  [G(u)-G*H(u)  ]e***^du 

Jo 


id^(t) 


m2(t) 


in(p^l)  ue*"^dG(u) 
f  [i-G(u)  ]e"“^u 

Jo 


^ot  at 
e  =  n^e 


m(iH-l)  I  ue*'^^dG(u) 


at  at 
e  =  n2e 


where  a  satisfies  /J”  e"^%G(u)  =  aud  so 
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Hg  = 


a[^~  m(p+l)] 


in(lH-l)  f 

Jo 


ue"°'^dG(u) 


5.9  Comparison  vlth  discrete  case. 

Checking  with  the  discrete  case  as  beforcj  ve  obtain 


M  = 


1 

|qd 


m(p+l)  J 


m“  = 


0 


[m(p+l)  f 


and 


„  qd 

mgTn)  m(p+l)-l 

for  both  the  discrete  matrix  case  and  the  "discretized"  continuous  case 
We  thus  obtain 
Lemma.  If 


(a^p)M  =  x^^(a,p)  =  Tn(p+i)(a,p)  , 

then 


qd  a 


00  « 


71 


5.10  E[N^(t)],  E[N2(t)]. 


Let 


B(s,t)  =  2  P[%(t)=J]s'^ 


J=0 


ft  n 

B(s,t)  =s[l-G(t)+/  [ph^(B(s,t-u))+(ih(B(s,t-u))JdG(u)] 

Jo 


E[No(t)]  ~ 


l_m(p+l)a  r  ue'^^dGCu)] 

Jo 


^0*  V,  cet 

e  =  . 


To  find  E[N^(t)],  set  H=0,  so  that  G^^Cs^^Sg, t)=s^.  With  this 
change. 


9G«(s  ,B„,t) 
E[W^(t)]=~2  1  2 


■Si: 


Si=82=1 


and 


E[N^(t)] 


“oi® 


at 


5.11  Convergence  in  qimdratlc  mean  and  asymptotic  correlations. 

The  arguments  are  as  in  the  previous  two  cell  type  model. 

Altered  steps  are  sketched-  As  before,  we  study  W^(t),  WQ^(t),  1—1,2, 
Define 

00 

B(s  ,S2,t,T)  =  2  P[Z,(t)=j;  Z,  (t+T)=k]s^sS 

j,k=0 
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/'t  2 

B(s^,S2,t,T)  =  1-G(t+T)+J  [ph  (BCs^^Sgjt-u^T)) 


+ 


qb(  1, t-u) )h (b( s^, Sg, t-u, T ) ]dG(u) 


4 


+  s^f  [ph  (GgCsgjljt+T-u))  + 


+  qb(G^(82flft+r~'u))h(G^(B^,l,t-i-r-u))]dG(u.) 


^  B(b- ^ 

E[Z^(t)Z^(t+T))  =  ^  ^ 


’1  '^®2 


where 


E[Z^{t)Z3^(t+T)]  ~  Ae^°*'*‘^ 


A  = 


[h"(l)  [l+p]+2pni^]  f  e"^'^^dG(u) 

Jo 


l-m(p+l) 


r 

Jo 


Define 


DO  ,  , 

C(B,,Sp,t,T)  =  2  P[Zp(t)=Ji  Zp(t+T)=k]s^s. 

^  i,k=0 

rt  2 

C(e^,Sp,t,T)  =  s^SgCl-GCt+t) ]  +J  [ph  (C(s^^S2,t-u,T)) 

+  qh(C(s^,B2>t-u,T)) ]dG(u) 
ft+T  2 

+  I  [ph  (GgCl, Sgjt+T-u))  +  qhCGgCljSg^t+T-u)) ]dG(u) 

•^t 


Tnfrr  {^\r7  / ,  .  _ M  2  .  2crt:+ccr 
ELZ2(i.;zj2\.u-r  t;  j  ~  ng  Ae 
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Define 


00 


D(s  ,S2,t,T)  =  X  P[N  (t)=j;  N.(t+T)=k]s^s; 
k>J>0-^ 

ft  2 

=  1-G(t+T)+J  [ph  (D(s^,S2,t-u,T)) 

+  qb(s^S2)h(DCs^^S2,t-u,T))]dGCu) 

rt+t 

+  /  [qD(s2)li(D(s2,l,t+T-u,0))+ph  (D(s2>ljt+T-u,0) ) ]dG(u) 


E[N^(t)N^(t+T)]  ~ 


Define 


H(s  ;,S2,tjT)  =  X  PtNp(t)=ji  Iirp(t+T)=k]s^Sp 

rt  2 

S2^tyT)=s^e2[l-G(t+T)+  j  [ph  (H(s^,S2,  t-u,T)  )+qh(H(sT ,  t-u,T) )  ]dG(u) 
rt+T 

+J  [ph  (H(s2^1,t+T-u,0))+qh(H(sg,l,t+T-u,0)) ]dG(u) 

E[Hg(t)N2(t+T)]  ~  nQ2 

These  results  suffice  to  show  that 

E[W^(t)-W^(t+T)]^  ->0  1=1,2 

E[WQ^(t)-WQ^(t+T)]^  -.0  1=1,2 

so  that  w^(t)  -»W^  and  1=1? 2, 


in  quadratic  mean. 


The  results  also  yield  that  the  limiting  pairwise  correlations 

among  W^(t),  W^Ct),  and  W^^Ct)  anproaeh  1  as  t-.»,  since 

we  have 

Lemma. 

A-1  =  Ih"(l)[l+p3  +  2pm^]c  Q 

l-m(p+l)c 

where 


Proof. 


c  >  [ 


Jo 


e-™dG(u)]2 


mtp+lj 


It  suffices  to  shov 


[h"(l)[l+p]  +  2pni^  +  m(p+l)][^^i^]^  -1  >  0 


or 


hut 


Since 

and 

it  follows  that 


h"(l)  +  m  >  m^(^) 
P+1 


h"(l)  +  m  >  m^  >  m^{^-i~)  . 


pH 


E[W^]  =  E[W  ]  =  1  1=1^2 


01 

E[W^]  =  E[Wq^]  =  A  i=l,2 


W 

1 


a.e. 
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APFEMDIX 


Lemma  1.  Consider  the  eq[iiatlon 


K(t)  =  f  (t)  +  m  [  K(t-u)dG(u) 
-*0 


G  is  a  distribution  function  on  (O,  «),  G(0)  =0,  G  is  not  a  lattice 
distribution,  f  is  a  known  function  bounded  on  every  finite  interval, 
m  >  1.  Define  a  >  0: 


[  e"^dG(u)  =  i  - 

Jo  ™ 

Suppose  f(t)e”^^-^0  as  t  and  f(t)e"^^  is  bounded  and 

integrable  in  (O,  »),  then 


Lemma  2.  In  the  equation 

K(t)  =  f(t)  +  m  r '^K(t-u)dG(u) 

h 

Suppose  m  <  1,  and  G  is  a  distribution  function,  G(o)  =0.  G  may 
or  may  not  be  a  lattice  distribution.  If 


then 


lim  f (t)  =  c, 
t-^  » 

K(t)  SB  t  -»oo  . 

^  1-m 


These  lemmas  are  from  [3],  Ch.  6,  appendix. 
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